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Abstract 

The Green function of the quark-antiquark system in the confin- 
ing background field is analysed using the Feynman-Schwinger for- 
malism. The Hamiltonian for the case of massive spinning quarks is 
obtained in the form containing essentially nonhermitian part. The 
eigenvalue problem for such type of the Hamiltonian is discussed, and 
it is shown that no complex eigenvalues arise. The corresponding 
nonunitary Foldy-Wouthuysen transformation is performed to obtain 
the hermitian Hamiltonian, and the standard spin-orbit interaction 
term is recovered. 
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1 Introduction 



The most natural way to discuss quark confinement at the constituent level 
is in terms of Wilson loop 

W{C) =Tr P exp ig j> A^X^dz^,. (1) 

The area law asymptotics for the Wilson loop (1) averaged over the vac- 
uum gluonic field 

< W{C) >~ exp{-(TS), (2) 

leads to the linear potential between heavy quarks, while the essentially non- 
local and velocity- and spin-dependent interaction generated by nonpertur- 
bative QCD exhibits itself as corrections of order of to the leading linear 
confinement potential term. These corrections, known as Eichten-Feinberg- 
Gromes relations were derived in several ways, see [1-3]. 

In the present paper we discuss the approach based on the most straight- 
forward use of the Feynman-Schwinger representation for the quark Green 
function [3-6]. The Green function of the quark in the given external field is 
the product of the quadratic Dirac propagator and the linear Dirac operator: 

—^ = {m + D) ^ (3) 

m — D m? — D'^ 

In the approach under discussion the Feynman-Schwinger representation 
is written out only for the quadratic part of the expression (3), and, as the 
result, the effective Hamiltonian of the qq system contains the nonhermitian 
spin-dependent part. The energy eigenvalues, however, are real. Our aim 
is to study this unusual phenomenon, and to demonstrate that the proce- 
dure is quite legitimate at least in leading order in l/w?, and the resulting 
interaction is equivalent to the standard one [1, 3]. 



2 qq Green function and the effective Lagrangian 

First we consider the Green function of quark- antiquark system [4]. 

Gqq^lv<x\ ^—^\y>V<y\ ^-^|a; > } (4) 

\ mi- Di 1712- D2 I A 
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where D is Dirac operator and averaging is performed over a nonperturbative 
field A. In what follows we assume m2 — > oo. For finite 1712 the result is 
generalized wihout difficulties. 

Each operator ^^^3- is represented in the form (m + .D) ^^^^^ , to expo- 
nentiate the positively defined value of the denominator( Feynman-Schwinger 
method [5]). For the latter we have 



< X 



(5/ 

where = d^ — igA^, P is an ordering operator and, in the Euclidean space, 

= ^ [t^x, 1v\ ; {in, Iv} = 25^^. (6) 
Linear operator D{A) corresponds to — |i^(s)7^ in the functional integral. 



[4, 7]. Indeed, the matrix element < x \ \ y > can be represented as 

00 

j ds < x \ -iDfj,e'^^' \ y > = 


00 

lds<x\{p,-gA,)e-^P-^^^'^\y> , 


where the linear operator —iD^ = P^l ~ qA^^ acts on the final state < x \; 
introducing the functional integration over the momenta one has 

7 r -jdrmT)-gA{z{T))f+ipz} 

j ds j DzDp{p^{s) - gA^{s))e 



JipzdTl § ^ - J(p-gAfdT-j gEFdT 



= JdsjDzDpeo 2^k^.e 

Integration by parts of this expression yields the replacement Df^ —\zn \r=s 
So we are spared from the necessity to average preexponent factor D{A) 
at the cost of appearance of in the preexponent in the functional integral. 
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Since this linear factor belongs to an ending point(x or y), we consider the 
exponent only. 

Integration over p gives the kinetic term + ^, Wilson loop W{C) — 
P exp ig J dz^A^ and the spin-dependent term, [4] . 

Instead of integrating over the closed contour we integrate over the area, 
using the cluster expansion [4]: 

{W{C))j^ = ^■^vY.^SwT j da^,^,...da^^,„ « F^,,,(l)...F^„,„(n) » 

(7) 

where da^i, = a^^dfl is the area element and we run over all the points with 
coordinates w^{[3,t) inside the contour Wfj,{l3,T) — Zn{t)(3 + Zij_{t){l — /?), 
< /9 < 1, = rdjddt] ra^u = w^w^ — w^w^. 

Since < T.F{x)W{C) >a= ^fiuj^^ < W{C) >, we have exponent 
exp{S^,^^^} which is really a shift operator: 

(Operator E^,^ is defined on the trajectory z (or z)). For bifocal correlators 
we have Wilson loop average expression: 

< W{C) >= exp{(-/y^ « Fi,^{w)F^p{w') » a^^(^/;)a«/3(w')dQdQ')}, 



and at r > rg,r > Tg < W{C) >= ex.p{-a J dtdpryfa?/2}, where a is 

oo 

defined from a — Cirg J D{x'^)xdx and the constant C and function D are 



taken from the expression for bifocal correlator [3] : 

<< F^^{w)Fccp{w') >>= const ([S^a^vp - 5^p5va)D{h?)+ 
+ l/2(9c,(/i^5;3^ - hjp^) + dpihJa^c - hf,5^^))Di{h^)) ,h^w-w' 

For spin-orbit term we obtain the expression 

(T Ea .0" an + av± _ , 

r (8) 



2/^ 2// ^ 



^2 
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in the Euclidean space (Instead of integration over zq we integrate over /i in 
the functional integral, [6]: 1/ {2ii{zo)) = ^,Vi = j^), fi plays the role of an 
effective mass. 

In the Minkovsky space (8) has the form 

a ian + av± 
and the effective Lagrangian equals to 



instead of tfj_ in 



In the case of finite mass m2 we get v^f^ (3 + vf^\l — (3) 

the area law, [8] , and should replace by for the spin-orbit terms. 

How do these expression change if high order correlators are taken into 
account in (7)? At distances r ^ Tg the same expressions are vahd and the 
only contribution from other correlators is reduced to the renormalization of 
the string tension a. Only Kronekker part of correlators( D for the bilocal 
one) gives the contribution to a since the parts with derivatives (like Di) 
are suppressed at large distances. For small distances r <^ Tg high order 
correlators give corrections {r/TgY. 

The spin-orbit expression also contains Tg/r terms( see for heavy masses 

[3]) 



3 Effective Hamiltonian 

The expression (10) for the Lagrangian makes sense only in the context of the 
Feynman-Schwinger representation, because the last term in (10) contains 7- 
matrices. To pass to Hamiltonian formulation one has to define the canonical 
momentum as p = Clearly, such procedure is possible only for heavy 
quark, when the last term in (10) is treated as perturbation; otherwise the 
expression for momentum would contain 7-matrices. 

For heavy quark the area-law term becomes ar J d(5\Jl — v\(5'^ ~ ar{l — 
\v']_), integration over 11 gives 11 = m, the Hamiltonian has a form 
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(11) 




(12) 



and we are left with the nonhermitian part in the Hamiltonian. The appear- 
ance of it is not surprising, it is caused by the fact that we have exponentiated 
only the quadratic part of the Dirac operator, and the projective operator 
has not been exponentiated. 

It appeares, however, that the Hamiltonian (11) has real eigenvalues. 
There is nothing mystical in such situation, and the corresponding examples 
are given in the Appendix. In the Hamiltonian (11) Vsl is considered as the 
perturbation. Taking the wave function of zero approximation as that of a 
free massive particle [££/), we get the matrix element of ian in the form 



It is necessary to work carefully with new wave functions if we want 
to apply perturbation theory for new Hamiltonian: new metric norm M is 
introduced: 



so that (HM)^ = HM (the pseudohermitian condition [9]) 

The Hamiltonians of such kind were considered in detail in the paper [9]. 
According to the theorem in [9] we can transform pseudohermitian Hamilto- 
nian with real eigenvalues to the hermitian one, and vice versa, if we obtain 
a hermitian Hamiltonian by a nonunitary transformation we can expect real 
eigenvalues of the original Hamiltonian. 

Using Foldy-Wouthuysen transformation, [8] 




2m 




mr 




ip — exp {iS)il), H — exp {iS)Hexp {—iS) 



with 




we get 



6 



H = m+^ + ar--^ + VsL (13) 

V — ^ _|_^|— (\^\ 

2mr 2mr j Am?r Am?r^ 

which is well-known, [1-3,8]. The contribution from the electric field reduces 
the one from the magnetic field by the factor 2. 




4 Discussion and conclusions 

The nonunitary Foldy-Wouthuysen transformation which leads to equations 
(13), (14) should be compared with the usual one. If one first applies the 
usual unitary Foldy-Wouthuysen transformation for the quark in the given 
external field, writes out the Fcynman-Schwinger representation for the trans- 
formed Green function and averages over the background field after that, as 
it was done in [3], no problems with nonhermitian Hamiltonian arises, and 
one arrives to the expressions (13), (14) straightforwardly. It should be no- 
ticed that the final Hamiltonian and the corresponding equations are not 
equivalent to Dirac equation for one particle in an external field. This is 
the case since after the averaging over nonperturbative field we get nonlo- 
cal interaction, and the real dynamical object is a string with quarks at the 
ends. The effective field is distributed between quarks and effective string. 
The spin-dependent interactions, on the contrary, are local and are defined 
along the contour: it is quark (antiquark) which has spin, and feels the dy- 
namics. We have got area law and spin-orbit term (14) using successively 
Feynman-Schwinger formalism without introduction of the effective field for 
a particle. On the other hand, if one wishes to consider the problem of a Dirac 
particle in an external field, proceeding from the intermediate result - area 
law, and, in the case of a heavy quark, from the specific form of the string 
correction, which can be treated as external potential A = l/3n x L, [10], 
doing in such a way one should obtain with neccessity 1/6 in (14) instead 
of 1/4, [11]. At this point the question arises: why this external potential is 
treated as something that should be substituted into the linear Dirac equa- 
tion? As the original Feynman-Schwinger representation was used only for 
the quadratic part of Dirac operator, this procedure is not well grounded at 
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all. 

We have demonstrated that it is possible to derive the spin-orbit force in 
the Wilson loop context, the results coincide with the well-known ones, and 
the Gromes relations are satisfied. We show, on the other hand, that if the 
particular form of Feynman-Schwinger representation is used, in which the 
projective operator is not inchidcd into the path integral, such procedure 
does not allow to go beyond the leading order in l/rn?. 
The auther is grateful for useful discussions to A.Yu. Dubin, Yu.S. Kalash- 
nikova, Yu.A. Simonov and acknowledges financial support of the Russian 
Fund for Fundamental Research N. 96-02-19184a. 



Appendix. Some examples of pseudohermitian 
Hamiltonian operator 

The first example is Klein-Gordon equation for a free heavy-mass particle, 

[9]: 

(Po -P^ - m^)(t> = 0, 
which can be transformed into: 

{^Jf^^^ + E){^Jf^^^ - E)(t) ^ Q 

and after expanding the root becomes 

2 2 4 

{m+ ^ + E){m+ ^ - E)(t> - -^0 = 
2m 2m Aw? 

2 

Defining % = — 0, we obtain the equation Hil) — Eijj for the 

two-component wave function ^jJ — with the following Hamiltonian: 

H has real eigenvalues although it is nonhermitian operator. For such Hamil- 
tonian another norm M should be defined: 
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^ 1, 

j ip+Mip = 1. 

The second example is Klein-Gordon equation for a massive particle in an 
external field: 

where Pq = -ido - Ao , Pi ^ Pi - Ai. 



For new wave function = (Pq + ^)(f> , ^ — + w? we have the 
equation: 

p,-n+[n,Po]p^^4> = Q (A.2) 

For heavy mass ^ ~ + ^ — 
and 



(A.3) 

The first order term is nonhermitian but after the nonunitary Foldy- 

Wouthuysen transformation = expiS(j), with 5" = ~'^P'^ ^^e new Hamil- 
tonian becomes a hermitian one up to the second order of 1/m: 

As we can see nonutary transformation was required in order to represent 
the Hamiltonian in a hermitian form. 

The third example is quadratic Dirac equation for a particle in the Coulomb 
field: 

- m^)^ = 0. 

(P2 -f-m"^- y.F)iIj = 0. (A.5) 
Since PqiP = {E — Aq)iP, {Aq = a/r,a < acru) we have the equation: 
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{m + e- Aq) ip = [p"^ + m"^ + SFJ ijj 
So the final equation is Hip = Eip, where 

H = m + Ao + — ^ ^ —. 

2m 2m 2m 



For a nonrelativistic particle it is possible to replace e — Aq hy ^,and 



H = ,n + A, + f-f,+ '^ (A.6) 

2m am-^ 2m 

The last term is nonhermitian although from the linear Dirac equation we 
have real eigenvalues for energy. 
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